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摘 要 : 叶 面 积 指数 (LAI) 是 描述 植被 生长 状况 和 冠 层 结构 的 一 个 重要 参数 ， 快 速 获取 大 面积 植被 与 作物 
LAI 对 于 生态 系统 科学 研究 、 农 林业 生 产 指导 具有 十 分 重要 的 理论 和 实践 意义 。 本 研究 选取 海南 岛 典 型 热带 
作物 一 一 橡胶 树 为 研究 对 象 ， 构 建 基于 卫星 遥感 植被 指数 的 橡胶 林 LAI 佑 算 模 型 并 分 析 其 变化 规律 。 结 果 表 
明 ， 相 较 于 归 一 化 植被 指数 (NDVI) 、 绿 色 归 一 化 植被 指数 (GNDVI)、 比 值 植被 指数 (RVI) 和 宽 动 态 范 
围 植 被 指数 (WDRVI) 四 个 指数 ， 增 强 植 被 指数 (EVI) 、 土 壤 调节 植被 指数 (SAVI) 、 差 值 植 被 指数 
(DVD 和 改良 土壤 调节 植被 指数 (MSAVI) 四 个 指数 同 LAI 之 间 的 相关 性 较 高 。 构 建 的 基于 不 同 植被 指数 的 
橡胶 林 LAI 估 算 模 型 (一 元 线性 、 指 数 和 对 数 模型 ) 中， 基于 EVI 指 数 的 橡胶 林 LAI 一 元 线性 估算 模型 效果 
最 佳 ， 其 决定 系数 已 为 0.69。 经 验证 ， 该 模型 估算 的 橡胶 林 LAI 精 度 较 高 ， 观 测 和 模拟 的 橡胶 林 LAI 线 性 拟 
合 尼 为 0.67， 均 方 根 误差 RMSE 为 0.16， 平 均 相 对 误差 RE 为 -0.25% ， 但 在 橡胶 林 LAI 中 值 区 存在 低估 现象 ， 
同时 在 LAI 高 值 区 和 低 值 区 存在 一 定 的 高 估 现 象 。 从 空间 分 布 来 看 ， 海 南 岛 橡胶 林 LAI 高 值 区 (4.40~6.23 ) 
主要 分 布 在 海南 岛 西部 做 州 、 白 沙 等 市 县 ，LAI 中 值 区 (3.80—4.40) 主要 分 布 在 海南 岛 中 部 澄 迈 、 证 昌 、 
琼 中 等 市 县 ，LAI 低 值 区 (2.69~3.80) 主要 分 布 在 海南 岛 东 部 和 南部 的 定安 、 琼 海 、 万 宁 、 乐 东 、 三 亚 等 
市 县 。 总 之 ， 构 建 的 基于 EVI 指 数 的 橡胶 林 LAI 一 元 线性 估算 模型 精度 较 高 ， 克 服 了 NDVI、GNDVI、RVI 等 
植被 指数 容易 出 现 指数 饱和 问题 ， 具 有 较 好 的 科学 性 和 良好 的 推广 应 用 价值 。 
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1 引言 它 与 植被 的 许多 生物 、 物 理 过 程 (如 光合 作用 、 
蒸腾 作用 等 ) 有 关 ， 也 是 陆地 表面 过 程 模 型 (如 

叶 面 积 指数 (Leaf Area Index, LAI) 是 描 作物 生长 模型 、 碳 循环 模型 等 ) 常用 的 输入 参量 

述 植被 生长 状况 和 冠 层 结构 的 一 个 重要 参数 趾 ， ”之 一 外， 影响 着 陆地 生态 系统 的 物质 和 能 量 循 环 
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过 程 。 快 速 准 确 地 获取 区 域 乃至 全 球 尺 度 植被 
LAI 对 于 生态 系统 研究 具有 重要 的 科学 意义 。 
LAI 也 是 农作物 、 林 木 和 果树 等 生理 生态 、 良 种 
选 育 、 作 物 长 势 、 作 物产 量 评估 和 病害 评价 研究 
的 一 个 重要 参数 5。LAI 通 常用 于 定量 表征 作 
物 叶 片 的 玻 密 程 度 和 冠 层 结构 特征 ， 作 物 生育 过 
程 中 的 LAI 序 列 能 反映 作物 长 势 动态 变化 ， 可 用 
于 作物 生长 诊断 和 提高 作物 估 产 精度 o LAIA 
是 林木 和 林 分 群体 结构 合理 性 的 重要 标志 之 一 ， 
合理 的 LAI 是 充分 利用 光 能 、 保 证 林 分 高 产 优 质 
的 主要 条 件 '。 林 分 LAI 研 究 可 以 为 合理 栽培 和 
丰产 林 培 育 提供 理论 依据 ， 并 成 为 衡量 林 分 质量 
的 重要 指标 中。 测定 叶 面 积 的 准确 与 否 ， 直 接 影 
响 到 对 林木 和 林 分 合理 受 光 群 体 结构 的 研究 和 光 
能 利用 率 及 净 同 化 率 的 计算 等 "。 快 速 获 得 大 面 
积 作 物 LAI 对 于 指导 农林 业 生 产 具 有 十 分 重要 的 
实践 意义 。 

目前 ，LAI 测 量 与 估算 方法 大 致 可 分 为 直接 
测量 法 和 间接 测量 法 两 类 。 直 接 测量 法 包括 格 点 
法 、 长 宽 校正 法 、 描 形 称 重 法 、 鲜 重 打 孔 法 、 标 
准 枝 法 、 面 积 收获 法 、 落 叶 收 集 法 和 面积 仪 测定 
法 等 ”。 间 接 测量 法 包括 点 接触 法 、 消 光 系 数 
法 、 经 验 公 式 法 、 光 学 仪器 法 和 遥感 方法 ， 
其 中 遥感 方法 包含 光学 模型 反 演 法 … 和 统计 模 
型 法 O 等 。 传 统 的 LAI 地 面 直 接 测 量 通常 存在 
对 作物 和 植被 具有 破坏 性 、 耗 时 耗 力 、 并 且 获 得 
言 息 有 限 、 难 以 扩展 到 大 区 域 ， 不 能 呈 面 状 分 布 
等 局 限 5。 与 传统 直接 测量 方法 和 光学 仪 吉 方 
法 比较 ， 卫 星 遥 感 数据 具有 覆盖 面积 大 、 更 新 周 
期 得、 花费 相 对 少 等 特点 和 无 损 、 快 速 、 大 面积 
应 用 等 优点 ， 故 卫星 遥感 方法 成 为 无 损 、 快 速 、 
大 面积 估算 植被 LAI 切 实 可 行 的 实用 方法 e 

当前 ， 国 内 外 关于 LAI 遥 感 估 算 研 究 主 要 集 
mE AO UT AH, PA BK RL) 
BE EX 等 大 宗 农作物 的 研究 ， 极 少 关 注 
热带 作物 ， 如 橡胶 林 、 油 棕 林 、 木 莫 、 热 带 果 树 
等 LAI 的 遥感 研究 ”YW。 虽 然 胡 炊 华 等 "研究 
了 海南 不 同 品系 胶 园 的 LAI 及 其 变异 ， 并 对 实生 


树胶 园 的 最 适 LAI 作 了 分 析 ; Rusli Ñ Majid |’ 
开展 了 小 流域 橡 股 林 和 油 棕 林 的 LAI 制 图 和 监测 
TE; Chen 等 "利用 卫星 数据 估算 了 海南 做 州 
试验 农场 橡胶 林 LAI， 但 是 这 些 研究 范围 均 只 涉 
及 较 小 空间 尺度 ， 尚 未 扩展 到 大 面积 橡胶 林 LAI 
遥感 估算 。 

橡胶 树 (Hevea brasiliensis) 是 天 然 橡 胶 生 
产 的 主要 来 源 ， 其 广泛 种 植 于 热带 地 区 。 橡 胶 种 
植 业 不 仅 是 海南 农业 的 支柱 产业 ， 也 是 农民 经 济 
收入 的 重要 来 源 。 开 展 典 型 热带 作物 一 一 橡胶 树 
LAI 遥 感 佑 算 对 海南 岛 橡胶 种 植 空间 优化 、 长 势 
监测 、 产 量 估算 也 有 重要 意义 。 因 此 ， 本 研究 选 
取 海 南 岛 橡胶 林 为 研究 对 象 ， 构 建 基 于 卫星 遥感 
数据 提取 的 各 类 植被 指数 数据 和 野外 观测 的 橡胶 
林 LAI 数 据 的 各 类 经 验 估算 模型 ， 并 评估 各 类 模 
型 的 表现 ， 选 出 最 优 估算 模型 ， 对 海南 咏 橡 胶 林 
LAI 进 行 遥 感 估 算 ， 以 期 为 橡胶 林 长 势 遥 感 监 
测 、 产 量 遥 感 评估 提供 科学 参考 。 


2 数据 与 方法 


2.1 研究 区 概况 


选择 海南 岛 为 研究 区 ， 橡 胶 林 分 布 及 LAI 观 
测 点 如 图 1 所 示 。 海 南 岛 地 处 18? 10'N—20?10'N 
和 108*37E~-111?03 中 之 间 ， 岛 屿 轮廓 形似 一 个 
椭圆 形 “ 大 雪梨 ”,， 长 轴 呈 东北 至 西南 向 长 约 
290 km， 西 北 至 东南 宽 约 180 km， 海 岸 线 总 长 
1944 km， 总 面积 为 33,900 km ， 是 中 国 仅 次 于 
台湾 岛 的 第 二 大 岛 。 海 南 岛 地 势 为 中 部 高 、 四 周 
低 ， 气 候 属 热带 季风 海洋 性 气候， 全 年 暧 热 ， 年 
平均 气温 高 ， 气 温 年 差 较 小 ， 雨 量 充沛 ， 干 湿 季 
AG, Sh. AL KERES, M E, RFA 
REKEM. MFP D K THE a 2 XU N/E 
林 ， 热 带动 植物 资源 十 分 丰富 。 海 南 是 中 国 重 
要 的 天 然 橡胶 树种 植 地 区 ， 天 然 橡 胶 是 海南 第 
一 大 经 济 作物 。 根 据 2020 年 海南 省 统计 年 鉴 数 
H, 截止 到 2019 年 底 , 海南 天 然 橡 胶 种 植 面积 
526,900 ha， 占 全 国 种 植 总 面积 的 46.00%; FFH 
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面积 381,300 haa， 占 全 国 开 市 总 面积 的 48.56%; 
干 胶 产 量 330, 800 t， 占 全 国 干 胶 总 产量 的 
40.8596; 干 胶 产 值 33.77 亿 元 ， 占 全 国 干 胶 产 值 
AY 41.22% P", 


2.2 野外 实测 数据 


本 研究 在 海南 岛 主 要 植 胶 地 区 选择 有 代表 性 
的 典型 橡胶 树林 分 样 地 LAI 观 测 点 77 个 〈 图 1)， 
采用 对 角 线 数据 采集 法 ， 以 及 手持 GPS、LAI- 
2200 实 地 观测 随机 样 方 LAI (图 2)， 从 而 获得 橡 
胶 林 观测 林 分 样 地 LAI 及 其 相关 数据 。 为 了 与 陆 
地 卫星 遥感 数据 空间 分 辨 率 保 持 一 致 ， 随 机 样 方 
的 大 小 设置 为 30 mX30 m; 同时 为 了 分 析 橡 胶 
林 生 长 最 佳 时 期 (8 一 9 月 份 ) 的 LAI， 本 研究 选 
择 橡 胶 林 LAI 观 测 日 期 为 2017 年 8 月 15 日 一 9 月 
15 日 。 此 外 观测 时 间 选 择 在 每 天 上 午 9 一 12 时 之 
间 ， 以 便 与 卫星 过 境 时 间 相 匹配 。 统 计 77 个 橡 
胶 林 LAI 观 测 点 数据 ，2017 年 8—9 H LAI 平均 
值 为 4.33+0.27， 最 大 值 为 5.12， 最 小 
值 为 3.68。 将 地 面 观测 LAI 样 本 数据 数 
量 随 机 分 成 7:3， 分别 用 于 LAI 佑 算 模 
型 的 构建 样本 (54 个 ) 和 验证 样本 
(231). 


2.3 海南 岛 橡胶 林 空 间 分 布 数 据 


本 研究 使 用 的 海南 岛 橡 胶 林 空间 分 
布 数据 来 源 于 中 国 热带 农业 科学 院 橡 
胶 研 究 所 (http: //www.catas.cn/xjs/) , 
该 数据 是 利用 日 本 JAXA 的 PALSAR/ 
PALSAR-2 和 Landsat 等 长 时 序列 遥感 
影像 ， 协 同 橡胶 林 物 候 特 征 SENSORS 
定 的 橡胶 林 识 别 算法 ， 获 取 2015 年 海 
南 岛 橡胶 林 空 间 分 布 信息 ， 总 体 分 类 精 
度 >90% 2 。 


30m 
2.4 遥感 数据 来 源 及 预 处理 (b) 随机 样 方 


Lass 
Em iH 


i: 中 国 地 图 基于 自然 资源 部 标准 地 图 服务 网 站 GS(2019) 
1697 号 标准 地 图 制作 、 海 南 岛 地 图 基于 GS(2019)3333 号 标准 
地 图 制作 , 底 图 边界 无 修改 

图 1 2015 年 海南 岛 橡 胶 林 分 布 和 LAI 观 测 点 示意 图 


Fig. 1 Spatial distribution and LAI observation points of rubber 


plantations in Hainan Island in 2015 


Co) 橡胶 林 观 测 点 (d) LAI 2200 


P FA 注 : 绿色 矩形 为 样 地 内 随机 样 方 , 浅 绿色 国 点 为 随机 样 方 数 据 采集 点 
本 研究 选取 覆盖 研究 区 范围 的 iki 图 (b) 中 ,绿色 撼 形 为 样 地 内 随机 样 方 , 浅 绿色 圆 点 为 随机 样 方 数 据 采 集 


2017 年 8 一 9 月 的 Landsat 8 OLI 时 间 序 


图 2 橡胶 林 野 外 观测 示意 图 


Fig. 2 Sketch map of rubber plantation LAI field observation 


列 卫 星 和 遥感 影像 23 景 ( 表 1)， 并 对 卫 
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星 遥 感 数据 进行 几何 校正 、 辐 射 校正 、 大 气 校正 
等 预 处 理 。 选 用 Landsat-8 卫星 的 2、3、4、5、6 
和 7 等 6 个 波段 参与 植被 指数 计算 ， 影 像 空间 分 
辩 率 为 30 m。 由 于 海南 岛 属 热带 季风 海 详 性 气 
候 ， 很 难 获得 完全 无 云 或 低 云 的 影像 。 为 解决 多 
云 履 盖 对 研究 的 影响 ， 利 用 像 元 级 云 量 提取 算法 
(具体 原理 参考 文献 [23] ) 计算 研究 时 段 内 每 
一 景 影像 的 云 量 ， 结 合 研究 区 域 特征 和 季节 特征 
综合 评定 每 一 景 影像 的 得 分 ， 生 成 研究 区 内 逐 像 
元 云 量 最 小 影像 〈 图 3)， 合 成 用 于 LAI 佑 算 的 遥 
感 影像 数据 ， 以 此 克服 研究 区 多 云 多 雨 气 候 特征 
对 LAI 佑 算 结果 的 影响 ， 提 高 LAI 估 算 的 可 
靠 性 。 


108°30E 


109E 


109%30E HOWE 


HOSE 


注 :影像 中 黑色 部 分 为 云 履 盖 去 除 处 理 之 后 的 结果 
图 3 基于 2017 年 Landsat-8 卫 星 数 据 的 研究 区 云 量 最 小 
合成 影像 
Fig. 3 Composited Landsat-8 image with minimum cloud 


cover in 2017 in study area 
2.5 遥感 植被 指数 计算 


基于 以 往 LAI 遥 感 估 算 研 究 中 向 用 的 植被 
指数 ， 本 研究 中 选取 了 8 种 卫星 遥感 植被 指数 
( 表 2)， 包 括 增强 植被 指数 (Enhanced Vegeta- 
tion Index, EVI) 后、 归 一 化 植被 指数 (Normal- 
ized Difference Vegetation Index, NDVI) ^, 2% 
色 归 一 化 植被 指数 (Green NDVI, GNDVI) ^*, 
比值 植被 指数 (Ratio Vegetation Index, RVI) 7", 


5&1 Landsat 卫星 遥感 数据 信息 


Table 1 Landsat remote sensing data information. 


编号 REG fr 日 期 时 间 
1 123 046 20174E8 H 11 H 10:59:07 
2 123 046 201749 H 12 H 10:59:12 
3 123 047 20174E8 H 11 H 10:59:31 
4 123 047 201749 H 12 H 10:59:36 
5 124 046 201748 H 02 H 1:05:14 
6 124 046 201748 A 18 H 1:05:20 
7 124 046 20174£9 H 03 H 1:05:22 
8 124 046 201749 A 19 H 1:05:25 
9 124 047 201748 H 02 H 1:05:38 
0 124 047 201748 H 18 H 1:05:44 
11 124 047 201749 H 03 H 1:05:46 
p 124 047 201749 H 19 H 1:05:49 
3 124 048 201748 H 02 H 1:06:02 
4 124 048 201749 H 03 H 1:06:10 
5 124 048 2017 年 9 月 19 日 1:06:13 
6 125 046 20174E8 H 09 H 11:11:28 
7 125 046 20174E8 H 25 H 11:11:32 
8 125 046 201749 H 10 H 11:11:33 
9 125 046 201749 H 26 H 11:11:39 
20 125 047 201748 H 09 H 11:11:52 
21 125 047 201748 H 25 H 11:11:56 
22 125 047 201749 H 10 H 11:11:57 
23 125 047 201749 H 26 H 11:12:03 


YE: 卫星 为 Landsat 8， 传 感 器 为 OLI 


土壤 调节 植被 指数 (Soi Adjusted Vegetation 
Index，SAVI) 的、 宽 动 态 范围 植被 指数 (Wide 
Dynamic Range Vegetation Index, WDRVI) 2 、 
差 值 植被 指数 (Difference Vegetation Index, 
DVI) ”和 改良 土壤 调节 植被 指数 (Modified 
Soil Adjusted Vegetation Index，MSAVI) ?", Jy 
减 小 空间 尺度 效应 ， 计 算 植 被 指数 时 利用 样 区 中 
心 点 及 其 周围 3X3 邻 域 像 元 的 地 表 反 射 率 的 平 
均值 。 

本 人 研究 中 对 Landsat 8 OLI 卫 星 有 适 感 数据 进行 
收集 、 预 处 理 、 最 小 云 量 合 成 和 植被 指数 计算 均 
是 利用 谷歌 地 球 纪 
https: //earthengine.google.com/) 进行 。 谷 歌 地 
ER S| SFE Google He [E BS] E A EK. BE B ER RI 
学 资料 (尤其 是 卫星 数据 ) 进行 在 线 可 视 化 计算 


敬 (Google Earth Engine, 


ChinaXiv 合 作 期 刊 


202302.00209v1 


chinaXiv 


和 分 析 处 理 的 云 计 算 平 台 ， 该 平台 能 够 存 取 卫 星 ”数据 资料 进行 处 理 分 析 ， 提 供 了 足够 的 运算 能 
图 像 和 其 他 地 球 观 测 me I 资料， 并 对 这 些 。 ”和 高 效 的 运算 效率 ”。 
表 2 植被 指数 及 其 计算 公式 
Table 2 Vegetation index and their formulas 
缩写 英文 名 称 中 文 名 称 计算 公式 
EVI Enhanced Vegetation Index 增强 植被 指数 241 EVI-2.5 x (pq, — Prea)! (pas, + 600,4 — T-5P bine + 1) (1 
NDVI Normalized Difference Vegetation Index 归 一 化 植被 指数 ”NDVIE( pj 7 Prea (Pir + Prea) 2 
GNDVI Green NDVI 绿色 归 一 化 植被 指数 2 GNDVIE(p = Doa (Prir + Pereen) 3 
RVI Ratio Vegetation Index 比值 植被 指数 '”7] RVI vip! P rea 4 
SAVI Soil Adjusted Vegetation Index 土壤 调节 植被 指数 23] SAVI=1.5 x (p, — Pred )/ (Pir + Prea + 0-5) 5 
WDRVI Wide Dynamic Range Vegetation Index ” 宽 动态 范围 植被 指数 [3 WDRVIR(0.2 x p, — Prog)! (0.2 X p, + Prea) 6 
DVI Difference Vegetation Index 差 值 植 被 指数 3 DVI — p= Pra 7 
MU ende 
MSAVI Modified Soil Adjusted Vegetation Index 改良 土壤 调节 植被 指数 5 7 
(2 x p, + 1) - 8*(p,, = Prea )/2 8) 
AUS 
2.6 LAI 估 算 模型 构建 与 验证 3 结果 与 分 析 


在 分 析 观 测 LAI 数 据 与 8 种 遥感 植被 指数 之 
间 的 相关 性 的 基础 上 ， 利 用 54 个 样 地 观测 LAT 
数据 与 8 种 遥感 植被 指数 ， 基 于 经 验 统计 方法 ， 
构建 基于 不 同 植被 指数 的 橡胶 林 LAI 佑 算 模型 
(一 元 线性 、 指 数 和 对 数 模 型 )， 并 对 模型 表现 进 
行 评估 ， 并 选 出 最 优 估 算 模型 。 不 同 模型 的 表现 
采用 决定 系数 (Coefficient of Determination, R?) 
进行 评估 和 优选 。 利 用 选 出 的 最 优 估算 模型 估算 
海南 岛 橡 胶 林 LAI， 获 得 大 面积 橡胶 林 LAI 的 空 
间 分 布 ， 并 利用 剩余 23 个 观测 LAI 数 据 对 估算 
结果 进行 精度 评价 和 验证 。 以 估算 值 与 观测 值 之 
间 线 性 拟 合 的 决定 系数 民 、 均 方 根 误差 (Root 
Mean Square Error, RMSE) 和 相对 误差 (Rela- 
tive Error, RE) (公式 (9) 和 公式 (10) ) 对 模 


型 的 表现 进行 验证 和 评价 。 
" (M.— Y 
muse. | 24 o) (9) 
N 
M.- O. 
RE = 一 - x 100% (10) 


其 中 ，YN 为 数据 个 数 ; M 和 0 分别 为 橡胶 林 
LAI 佑 算 值 和 观测 值 。 


3.1 植被 指数 与 LAI 的 相关 性 


为 了 更 好 地 进行 估算 模型 构建 ， 计 算 植 被 指 
数 与 观测 LAI 之 间 的 相关 系数 。 由 表 3 可 知 ， 
EVI、SAVI、DVI、MSAVI 同 LAI 之 间 的 相关 性 
较 高 ， dde de AME 0.81, 0.82 fil 
0.81， 均 通过 p<0.001 显著 性 水 平 检验 。LAI 与 
NDVI、GNDVI、RVI、WDRVI 之 间 的 相关 性 较 
低 ， 其 相关 系数 分 别 为 0.61、0.65、0.67 和 0.65， 
均 通 过 p<0.01 显著 性 水 平 检验 。 

表 3 橡胶 林 8 种 植被 指数 与 观测 LAI 相关 系数 


Table 3 Correlation coefficients between 8 vegetation 


indices and LAI of rubber plantation 


EVI NDVI GNDVI RVI SAVI WDRVI DVI MSAVI 


r 0.83" 0.61" 0.65! 0.67 0.81" 0.65” 0.82" 0.81" 
注 : * 表 示 通 过 p<0.01 显著 性 水 平 ，** 表 示 通 过 p<0.001 显著 性 


3.2 估算 模型 构建 与 优选 


通过 构建 基于 不 同 植被 指数 的 橡胶 林 LAI 佑 
算 模 型 (一 元 线性 、 指 数 和 对 数 模 型 )( 表 4)， 
在 一 元 线性 模型 中 ， 基 于 EVI 指 数 的 橡胶 林 LAI 
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估算 模型 效果 最 佳 ( 公 式 (040) )， 其 决定 系数 
R30.69, ts F3ETNDVI, GNDVI, RVI, SA- 
VI、WDRVI、DVI 和 MSAVI 指 数 的 估算 模型 。 
在 指数 模型 中 ， 基 于 SAVI 指 数 的 橡胶 林 LAI 估 
算 模型 效果 最 佳 (公式 (24) )， 其 决定 系数 尼 
为 0.67， 高 于 基于 其 他 植被 指数 的 估算 模型 。 在 
对 数 模型 中 ， 基 于 DVI 指数 的 橡胶 林 LAI 佑 算 模 
型 效果 最 佳 (公式 (31) )， 其 决定 系数 尺 为 
0.65， 高 于 基于 其 他 植被 指数 的 估算 模型 。 综 合 
比较 ， 选 择 基于 EVI 指 数 的 一 元 线性 模型 作为 橡 
胶 林 LAI 佑 算 的 最 优 模型 ， 并 估算 海南 岛 橡胶 
林 LAI。 
£4 基于 卫星 遥感 植被 指数 的 LAI 估算 模型 


Table 4 LAI estimation models based on satellite remote 


sensing vegetation indices 


植被 指数 估算 模型 R NH 
y=2.2582x+2.8622 (11) 0.69 

EVI — y-3.0762e5»* (12) 0.66 54 
y=1.3527In(x)+4.9312 (13) 0.63 
y=2.5583x+2.2212 (14) 0.34 

NDVI | y-2.6164e 6056: (15) 0.36 54 
y71.7281In(x)-4.6736 (16) 0.1 
y=4.1216x+1.2552 17) 0.44 

GNDVI _ y=2.0872e°75* (18) 0.46 54 
y-2.76511n(x)*5.1474 19) 042 
y=0.0645x+3.5703 20) 0.41 

RVI y-73.6154e9051s (21) 042 54 
y=0.5797In (x) +2.9245 (22) 0.38 
y=3.2588x+2.4713 23) 0.66 

SAVI — y-2.8047e" 759 (24) 0.67 54 
y=1.6999In(x)+5.2991 25) 0.63 
y=1.2421x+3.8685 (26) 0.38 

WDRVI  y=3.8740e°2°40 (27) 0.40 54 
y=0.1783In(x)+4.5408 (28) 0.15 
y=0.0004x+3.0470 (29) 0.65 

DVI — y-3.172e99008: (30) 0.65 54 
y-1.26971n(x)-6.0374 (31) 0.65 
y74.9797x42.3614 (32) 0.66 

MSAVI  y-2.7360e 56 (33) 0.67 54 

( 


y=1.8432In(x)+6.0540 


3.3 估算 结果 与 分 析 


利用 选 出 的 最 优 佑 算 模 型 估算 海南 岛 橡胶 
PK LAI， 获 得 大 面积 橡胶 林 LAI 的 空间 分 布 


(图 4)。 经 统计 ， 佑 算 的 海南 岛 橡胶 林 LAI 平 均 
值 为 4.32+0.18， 最 大 值 为 6.23， 最 小 值 为 2.69。 
从 空间 分 布 来 看 ， 海 南 岛 橡胶 林 LAI 高 值 区 
(4.40—6.23) 主要 分 布 在 海南 岛 西 部 人 黎 州 、 白 沙 
等 市 县 ，LAI 中 值 区 (3.80—4.40) 主要 分 布 在 
海南 岛 中 部 港 迈 、 屯 昌 、 琼 中 等 市 县 ，LAI 低 值 
区 (2.69—3.80) 主要 分 布 在 海南 岛 东 部 和 南部 
的 定安 、 琼 海 、 万 宁 、 乐 东 、 三 亚 等 市 有 具 。 海 南 
岛 橡 胶 林 LAI 空 间 分 布 存在 异 质 性 的 原因 包括 以 
下 几 个 方面 : 一 是 由 于 海南 岛 橡胶 林 大 部 分 年 龄 
低 于 40 年 中 ， 这 一 年 龄 段 的 橡胶 树 正 处 于 生长 
旺盛 期 ， 生 物 量 大 、 生 产 力 高 ，LAI 整 体 也 较 
高 ; 二 是 人 研究 区 年 龄 小 于 25 年 的 橡胶 林 大 多 分 
布 在 海南 岛 西 部 稀 州 、 白 沙 等 市 县 OU, GERE 
南 岛 橡 胶 林 LAI 高 值 区 域 主要 分 布 在 海南 岛 的 西 
部 市 县 ; 三 是 海南 岛 东 部 和 南部 易 受 台风 灾害 影 
响 ， 造 成 橡胶 林 损 害 ， 导 致 该 区 域 橡胶 林 LAI 值 
较 低 。 此 外 海南 岛 橡胶 林 主 要 种 植 在 海拔 100— 
600 m 山坡 地， 使 得 其 LAI 空 间 分 布 差 异 受 地 形 
影响 较 大 ， 分 布 于 山坡 地 区 不 同 海拔 、 坡 度 、 坡 
向 的 橡胶 林 接 收 到 的 太阳 辐射 呈现 空间 差异 ， 导 
致 植被 生长 受到 影响 ， 这 也 使 得 海南 岛 橡胶 林 
LAI 呈 现 空间 差异 。 


3.4 精度 验证 与 分 析 


比较 估算 和 观测 的 橡胶 林 LAI 结 果 (图 5) 
可 知 ， 橡 胶 林 LAI 佑 算 模型 精度 较 高 ， 模 型 估算 
结果 基本 处 于 95% 置 信 区 间 ， 观 测 和 模拟 的 橡胶 
Fk LAI 线 性 拟 合 决定 系数 尼 为 0.67， 通 过 p< 
0.001 显著 性 水 平 检验 ， 均 方 根 误差 RMSE 为 
0.16, 平均 相对 误差 RE 为 -0.25%。 但 在 橡胶 林 
LAI 中 值 区 (3.80—4.40) 存在 低估 现象 ， 同 时 
在 LAI 高 值 区 (>4.40) 和 低 值 区 (<3.80) 存在 
一 定 的 高 估 现 象 ， 这 可 能 是 由 于 海南 岛 橡胶 林 总 
体 处 于 生长 旺盛 期 〈 年 龄 小 于 25 年 的 面积 比例 
为 52.54% 7 ) ， 其 植被 指数 和 LAI 整 体 处 于 较 高 
值 ， 造 成 其 区 分 度 较 低 ， 使 得 橡胶 林 LAI 估 算出 
现 偏差 。 
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m 克服 了 研究 区 多 云 多 雨 气候 特征 对 LAI 估 算 结 
的 影响 ， 提 高 了 LAI 估 算 的 可 靠 性 。 基 于 到 感 大 
数据 平台 谷歌 地 球 引擎 构建 的 适 感 数据 集 ”， 
为 LAI 估 算 研 究 与 应 用 提供 了 充足 的 数据 基础 。 

橡胶 林 是 典型 的 热带 作物 ， 其 树 形 高 大 、 杆 
被 茂密 、 生 产 力 高 。 而 在 植被 茂密 、 生 产 力 高 的 
T 热带 地 区 遥感 研究 中 ,单一 植被 指数 往往 无 法 准 
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图 4 基于 EVI 指 数 和 一 元 线性 模型 的 研究 区 橡胶 林 LAI 
空间 分 布 图 
Fig. 4 Spatial distribution of rubber plantation LAI in the 


study area based on EVI index and linear model 


5 
R=0.67, RMSE —0.16 

MEAN, = 4.28 

45+ N=23 . 


估算 LAI 


45 5 


3 3.5 


4 
观测 LAI 
iki 灰色 阴影 部 分 为 95% 置 信 区 间 

图 5 研究 区 橡胶 林 LAI 估 算 结 果 验 证 


Fig. 5 Validation results of estimate rubber plantation LAI 


in the study area 
4 讨论 与 结论 


4.1 讨论 


海南 岛 地 处 热带 ,该 区 域 多 云 多 雨 ， 禾 盖 研 
究 区 的 光学 影像 数据 质量 往往 不 佳 ， 不 能 满足 研 
究 需 求 。 本 人 研究 中 ， 针 对 海南 多 云 多 雨 的 热带 气 
候 特 征 ， 采 用 像 元 级 最 小 云 量 影像 合成 方法 ， 获 
JX Landsat 时 间 序 列 不 同 波段 合成 的 纯净 影像 ， 


非 线 性 拉 伸 的 方式 增强 了 近 红 外 和 红 光 波段 反射 
率 的 对 比 度 ， 导 致 其 对 高 植被 区 具有 较 低 的 灵敏 
度 。 本 研究 中 ， 基 于 EVI 指 数 构建 橡胶 林 LAI 佑 
算 模 型 ， 较 好 地 解决 了 指数 饱和 问题 ， 这 是 由 于 
EVI 指 数 中 红 光 和 近 红 外 波段 的 范围 设置 更 罕 ， 
不 仅 提高 了 对 稀 玖 植被 探测 的 能 力 ， 而 且 减 少 了 
水 汽 的 影响 ， 同 时 ， 引 入 了 蓝光 波段 ， 对 大 气 气 
溶胶 的 散射 和 土壤 背景 进行 了 矫正 ， 可 以 稳定 地 
反映 热带 地 区 植被 的 情况 。 经 检验 ， 构 建 的 基于 
EVI 指 数 的 橡胶 林 LAI 佑 算 模 型 精度 也 较 高 ， 具 
有 较 好 的 科学 基础 和 实用 性 。 但 在 模型 推广 应 用 
过 程 中 ， 一 方面 需要 更 大 空间 覆盖 范围 的 地 面 样 
点 验证 ， 另 一 方面 还 需要 进行 长 时 间 序 列 的 地 面 
样 点 观测 ， 以 满足 未 来 LAI 时 空 反 演 和 时 空 分 析 
需求 。 

LAI 遥 感 估 算 方 法 主要 包含 光学 模型 反 演 和 
统计 模型 2 种 。 光 学 模型 反 演 法 是 基于 植被 的 双 
向 反射 率 分 布 函数 ， 是 一 种 建立 在 辐射 传输 模型 
基础 上 的 模型 ，LAI 作 为 输入 变量 ,采用 迭代 的 
方法 来 反 演 LAI。 这 种 方法 的 优点 是 有 物理 模型 
基础 ， 不 受 植被 类 型 的 影响 ， 然 而 由 于 有 些 模型 
过 于 复杂 ， 反 演 非 常 耗 时 ， 且 反 演 估算 LAI 过 程 
中 有 些 函 数 并 不 总 是 收敛 的 ， 从 而 限制 了 光学 模 
型 反 演 方法 的 使 用 和 发 展 "。 统 计 模 型 法 主要 
是 从 遥感 图 像 数 据 中 提取 植被 指数 与 野外 实测 
LAI 建 立 关 系 模型 ， 进 而 估算 植 被 LAI。 和 常用 的 
植被 指数 有 EVI 和 NDVI 等 ， 这 种 方法 在 不 同 植 
被 类 型 ， 如 森林 、 草 地 等 的 LAI 与 植被 指数 的 函 
数 关系 不 同 ， 在 使 用 时 需要 重新 计算 WMA 
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由 于 该 模型 方法 参数 输入 单一 ， 计 算 简 单 “” ， 
因此 ， 本 研究 也 采用 此 方法 进行 橡胶 林 LAI 的 遥 
感 佑 算 ， 也 取得 了 较 好 的 效果 。 

综 上 ， 未 来 需要 充分 利用 遥感 云 计 算 平 台 计 
ROE, FG UE, FEAT AEA Te 
源 的 组 合 优势 ， 发 展 更 高 精度 的 估算 模型 ， 如 辐 
射 传输 模型 、 机 理 模型 、 深 度 学 习 模 型 等 ， 提 高 
橡胶 林 等 热带 作物 LAI 估 算 精 度 ， 为 橡胶 林 长 势 
冯 感 监测 、 产 量 禹 感 评估 提供 科学 参考 。 


4.2 结论 


本 研究 通过 构建 基于 卫星 遥感 植被 指数 的 橡 
胶 林 LAI 估 算 模 型 ， 实 现 了 海南 岛 橡 胶 林 LAI 遥 
感 估 算 。 结 果 表 明 ， 相 较 于 NDVI、GNDVI、 
RVI 和 WDRVI 指 数 ，EVI、SAVI、DVI 和 MSA- 
VI 指数 同 LAI 之 间 的 相关 性 较 高 。 构 建 的 基于 不 
同 植被 指数 的 橡胶 林 LAI 估 算 模 型 (一 元 线性 、 
指数 和 对 数 模 型 ) 中 ， 基 于 EVI 指 数 的 橡胶 林 
LAI 一 元 线性 估算 模型 效果 最 佳 (公式 OD )， 
其 决定 系数 RR 为 0.69。 经 验证 ， 该 模型 估算 的 橡 
胶 林 LAI 精 度 较 高 ， 观 测 和 模拟 的 橡胶 林 LAI 线 
PEW RH 0.67, RMSE 为 0.16，R 为 -0.25%， 
但 在 橡胶 林 LAI 中 值 区 存在 低估 现象 ， 同 时 在 
LAI 高 值 区 和 低 值 区 存在 一 定 的 高 佑 现象 。 从 空 
间 分 布 来 看 ， 海 南 岛 橡胶 林 LAI 高 值 区 (4.40~ 
623) 主要 分 布 在 海南 岛 西部 人 黎 州 、 白 沙 等 市 
县 ，LAI 中 值 区 (3.80—4.40) 主要 分 布 在 海南 
Mm PRB. HAL Pea, LARE K 
(2.69—3.80) 主要 分 布 在 海南 岛 东 部 和 南部 的 定 
T., WE, AT. RAR, SWEAR. AO 
建 的 基于 EVI 指 数 的 橡胶 林 LAI 一 元 线性 估算 模 
型 精度 较 高 ， 克 服 了 NDVI、GNDVI、RVI 等 植 
被 指数 容易 出 现 指数 饱和 问题 ， 具 有 较 好 的 科学 
性 和 良好 的 推广 应 用 价值 。 
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Comparison of Remote Sensing Estimation Models for 
Leaf Area Index of Rubber Plantation in Hainan Island 


DAI Shengpei’”, LUO Hongxia", ZHENG Qian", HU Yingying", LI Hailiang", 
LI Maofen", YU Xuan", CHEN Banggian? 


(1. Key Laboratory of Agricultural Remote Sensing, Ministry of Agriculture and Rural Affairs, Beijing 100081, 
China; 2. Institute of Scientific and Technical Information, Chinese Academy of Tropical Agricultural Sciences/Key 
Laboratory of Applied Research on Tropical Crop Information Technology of Hainan Province, Haikou 571101, 
China;3. Rubber Research Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China) 


Abstract: Leaf area index (LAI) is an important index to describe the growth status and canopy structure of vegetation, is of 
great theoretical and practical significance to quickly obtain LAI of large area vegetation and crops for ecosystem science re- 
search and agricultural & forestry production guidance. In this study, the typical tropical crop rubber tree in Hainan Island was 
selected as the research area, the LAI estimation model of rubber plantation based on satellite remote sensing vegetation indices 
was constructed, and its spatiotemporal variation was analyzed. The results showed that, compared with correlations between 
LAI and the indices of normalized difference vegetation index (NDVI), green NDVI (GNDVI), ratio vegetation index (RVI) and 
wide dynamic range vegetation index (WDRVI), correlations were higher between LAI and the indices of enhanced vegetation 
index (EVI), soil adjusted vegetation index (SAVI), difference vegetation index (DVI) and modified soil adjusted vegetation in- 
dex (MSAVI). Among the LAI estimation models based on different vegetation indices (linear, exponential and logarithmic 
models), the linear estimation model based on EVI index was the best, and its coefficient of determination (R°) was 0.69. The ac- 
curacy of LAI estimation model was high. The linear fitting R^ of observed and simulated LAI was 0.67, the root mean square 
error (RMSE) was 0.16, and the average relative error (RE) was -0.25%. However, there was underestimation in the middle val- 
ue and overestimation in the high and low value area of LAI. The high LAI values (4.40— 6.23) were mainly distributed in 
Danzhou and Baisha in the west of Hainan Island, the middle LAI values (3.80 —4.40) were mainly distributed in Chengmai, 
Tunchang and Qiongzhong in the middle of Hainan Island, and the low LAI values (2.69 —3.80) were mainly distributed in 
Ding'an, Qionghai, Wanning, Ledong and Sanya in the east and south of Hainan Island. In summary, the linear estimation model 
for rubber plantation LAI based on EVI index obtained high accuracy, and has good values of popularization and appliance. 


Key words: leaf area index (LAI); rubber plantation; remote sensing; estimation model; Hainan Island 
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